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Using the time-dependent non-crossing approximation, we calculate the transient response of 
the current through a quantum dot subject to a finite bias when the dot level is moved suddenly 
into a regime where the Kondo effect is present. After an initial small but rapid response, the 
time-dependent conductance is a universal function of the temperature, bias, and inverse time, all 
expressed in units of the Kondo temperature. Two timescales emerge: the first is the time to reach 
a quasi-metastable point where the Kondo resonance is formed as a broad structure of half-width of 
the order of the bias; the second is the longer time required for the narrower split peak structure to 
emerge from the previous structure and to become fully formed. The first time can be measured by 
the gross rise time of the conductance, which does not substantially change later while the split peaks 
are forming. The second time characterizes the decay rate of the small split Kondo peak (SKP) 
oscillations in the conductance, which may provide a method of experimental access to it. This 
latter timescale is accessible via linear response from the steady state and appears to be related to 
the scale identified in that manner [A. Rosch, J. Kroha, and P. Woffle, Phys. Rev. Lett. 87, 156802 
(2001)]. 
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I. INTRODUCTION 

The theoretical predictions^^ of consequences of the 
Kondo effect for the steady state conduction through 
quantum dots began a decade ago. At low tempera- 
tures, a narrow resonance in the dot density of states 
can form at the Fermi level, leading to a large enhance- 
ment of the dot's conductance, which is strongly depen- 
dent on temperature, bias, and magnetic field. Many 
of these effects have been recently observed by a set of 
beautiful experiments by several groups4i£*£ These suc- 
cesses, supplemented by the anticipation that time de- 
pendent experiments are not far behind, have spurred a 
number of theoretical groups .LgaLiSiiiii^i to consider the 
effects expected when sinusoidal biases or gate potentials 
are applied. Indeed recent experiments^ have now seen 
Kondo sidebands. Also the predictions 14 of split Kondo 
conductance peaks have been observed in double 1 ^ and 
multiple 16 dots. Surprisingly, the application of steps or 
pulses, which can provide a less ambiguous measure of 
time scales than ac modulation, have been considered 
less extensively theoreticallyjiSiii 1 ^ and not at all to our 
knowledge experimentally. When pulsed voltage is ap- 
plied to the to the quantum dot level so that it suddenly 
is shifted into the Kondo regime, the conductance of the 
dot will begin to increase. The current saturates when 
the system reaches its new equilibrium configuration. In 
a previous investigation^ we considered a quantum dot 
biased by a small voltage and calculated the time de- 
pendent change in linear response conductance when a 
stepped potential was applied to a gate, thereby shifting 
the dot into the Kondo regime. Some general qualitative 



observations were made, which now can be made quan- 
titative through the study of a different configuration of 
voltage switching. 

In the present work we consider the response of a quan- 
tum dot, operating as a single-electron transistor, already 
subject to a finite dc bias when the dot level is shifted into 
the Kondo regime. When a finite bias is present across 
the dot, Kondo resonances appear at each of the leads. 
Due to the finite bias, these resonances are broadened 
compared to the zero bias situation. 19 While the equilib- 
rium or zero bias situation is relatively well understood?^ 
much less is know about the fully time-dependent situ- 
ation, with Refs. |2l| and 22 being most relevant to the 
present work. 

The present calculation determines within the non- 
crossing approximation (NCA), the transient current af- 
ter a gate pulse moves the level of the single electron 
transistor into the Kondo regime, under a large variety 
of temperatures and biases. We are able to extract two 
different timescales from our analysis, in addition to the 
trivial very rapid scale of the very small time response set 
by the tunneling rate to the leads (or the width of the 
virtual- level resonance). The first of these is character- 
ized by the time to reach a quasi-metastable point where 
the time-dependent conductance has essentially risen to 
its equilibrium value, and where the Kondo resonance has 
formed into a broad quasi-smooth structure of half- width 
roughly equal to the bias. The second is the timescale for 
the emergence and formation of the narrower individual 
split Kondo peaks (SKP), at a much slower timescale. 
This scale is also relevant for the decay of SKP conduc- 
tance oscillations, which may provide an experimental 
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access to it. 



II. MODEL 

We model the quantum dot by a single spin degener- 
ate level of energy edot coupled to leads through tunnel 
barriers. The Coulomb charging energy U prevents the 
level from being doubly occupied. 

The system may be described by the following Ander- 
son hamiltonian: 

H(t) = ^2e dot (t)n a + ^2 \tk n ka + {V k c\ a c a + H.c.) , 

a ktj 

with the constraint that the occupation of the dot cannot 
exceed one electron. Here c\ creates an electron of spin a 
in the quantum dot, with n a the corresponding number 
operator; c\ a creates an electron in the leads. 

For zero bias across the dot, the general features of the 
static equilibrium spectral density when the dot level e is 
sufficiently below the Fermi level are well known. There is 
a broad resonance of half-width T(e) = 2ir J^k | Vfc | 2 <5(e — 
efe) at an energy ~ edot- The notation T with no energy 
specified will always refer the value at the Fermi level. In 
addition, there is a sharp temperature sensitive resonance 
at the Fermi level (the Kondo peak), characterized by the 
low energy scaled Tk (the Kondo temperature), 

Tk = d (0^(-^), (2) 

where D is a high energy cutoff equal to half bandwidth 
when modeled by a symmetric flat band. Our calcula- 
tions here use a symmetric parabolic band of half band- 
width Do = 9r. We use D ~ Do/y/e, the choice that 
gives the correct normalization for the leading logarith- 
mic corrections in the Kondo model^ For our case where 
edot is in the band, E = D in Eq. J2J. Only for edot suf- 
ficiently below the band cutoff does the form E cc |edot|) 
expected from the Schrieffer- Wolff transformation from 
the Kondo model, result. For finite bias V, the Kondo 
peak splits into two sub-peaks^ at ±V/2 relative to the 
Fermi level, which we will always take to be at zero en- 
ergy. The nature of this splitting at large V has recently 
been elucidatediSi 

Our calculations use the non-crossing approximation 
(NCA), which is reliable for temperatures down to 
T < Tk 20 . The details of the time-dependent method 
of solution have been described in several previous 
publicationsi 2 ^ 2 ^ Throughout this work energies, tem- 
peratures, and biases are given in units of T, and times in 
units of 1/r, with h = fcs = e = 1. In the regions of pa- 
rameter space where Tk, T, and V are much smaller than 
edot and T, physical properties are functions of T/Tk, 
V/Tk and Xk£, alone, provided that the quantity under 
consideration is made dimcnsionlcss by appropriate fac- 
tors of Tk and Go = 2(e 2 /2-kK). In what follows we 



mostly focus on properties relevant to this universal re- 
gion, except when a property is directly relevant for ex- 
periments. 

In the present calculations, we investigate the transient 
electric currents trough the dot. The current into the dot 
depends on the time t as 

Jin(t) = ie^V k (cl a (t)c a (t)) +c.c. (3) 

ka 

It may be divided into contributions 7i e ft (t) and /right (i) 
by respectively restricting the k summations to the ap- 
propriate lead. For simplicity, we will only consider dots 
with left-right symmetry. The transport current is then 
= Tj[I\c{t(t) — /right (t)]- The finite bias on the leads 
is taken into account by introducing a time-dependent 
phase in V k in QJ. We calculate the Keldysh propaga- 
tors corresponding to the angular-bracketed expectation 
values in © for each lead, and hence obtain lit). 

III. TIME-DEPENDENT CONDUCTANCE 
RESULTS 

In this work our studies all start from an equilibrium 
steady state with a bias across the dot equal to V, but 
where the virtual level parameter edot (t) is at a negative 
value of magnitude sufficient that the initial Kondo tem- 
perature is much smaller than the physical temperature 
T, and the initial conductance is so small as to be negli- 
gible. In practice we use a starting value for time t < 
of — 5r. At t = this virtual level parameter is suddenly 
shifted to its final value, which we will simply call edot 
(with no time argument), with the bias V and virtual- 
level width parameter T unchanged. The majority of our 
calculations take edot to be — 2r, and we will denote the 
system so described as system one (SI). Its Kondo tem- 
perature Tk = 0.0022r. We also make a number of cal- 
culations where e do t = -2.225I 1 , for which T K = O.OOlir. 
This latter system we call system two (S2). 

When the dot level is in its lower position, the Kondo 
temperature of the dot is much smaller than the sys- 
tem temperature and the Kondo resonance is essen- 
tially absent. The spectral function is dominated by the 
broad virtual level of width ~ 2r centered roughly at 
e dot(* < 0). When the level is moved, a new virtual-level 
resonance of width ~ 2T is formed around the new dot 
level. The time scale for the formation of this resonance is 
riS The the Kondo resonance take longer time to form. 
In the G(t) curves shown in this section, only this Kondo 
time scale is apparent. 

A. Results for small bias 

We display first the results appropriate to very small 
V, i.e. V <C Tk, such that the conductance is in the 
linear response regime, for a variety of temperatures in 
system one. These are shown in Fig. ^ It is obvious 
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by inspection of the figure that the rise time of the con- 
ductance, however that should be defined, is increasing 
as the temperature is lowered. This increase saturates 
when the temperature gets to around Tk or lower. More 
quantitatively the inverse rise time (rise rate) is around 
7k at low temperatures, but increases with T for higher 
ones. These rates are much slower than the rate Y which 
is associated with the virtual level. 
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FIG. 1: Time dependent conductance vs. time for sys- 
tem one. The ordinate is the time dependent conductance 
G(t) = I(t)/V in units of the open-channel conductance 
Go = 2(e 2 /2-kTi). As before, I(t) is the time dependent cur- 
rent, and V the bias. The abscissa is the time t after the 
gate switches the dot into the Kondo regime, in units of 
The numbers in the legend are the temperatures T for the 
corrsponding curves in units of F. All curves are in the linear 
response region of approaching zero bias V. 



In Fig. [3 we combine results from system one and sys- 
tem two to test universality. We find the results sat- 
isfying. However, we believe that the small differences 
between the two sets of curves are not due to compu- 
tational inaccuracy, but rather are true deviations from 
universality, however tiny. The small time region Yt ~ 1 
is always non-universal, and around 25% smaller for sys- 
tem two than for system one. With our present algo- 
rithms, we cannot further separate the time scales using 
much smaller Kondo temperatures, so the results even 
on the Kondo time scale are still slightly affected by the 
nonuniversal fast contribution. For this reason, the con- 
ductances for system two will be a little smaller than for 
system one. In Fourier space, one could say that there 
is a small non-universal background, whose changes from 
system one to two are reflected in the final value of G. 
However, our results clearly show that the parts expected 
to be universal behave in this manner. 
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FIG. 2: Universality of time-dependent conductance curves. 
The SI curves are identical to the corresponding ones in Fig.|S| 
with the time axis rescaled to inverse Tk units, where here 
Tk = 0.0022r. The S2 are curves for system two (see text), 
which has a smaller Tk. 



B. Results for finite values of bias 

The finite bias calculations were performed in the same 
way, except that a constant bias V is present at all times. 
In Fig. we show the results for system one for a variety 
of different biases and two different temperatures. The 
simple descriptive facts that are evident by inspection 
of the rise-times implied by these curves are that: (i) 
for V <~ T there is little change; the rising rate propor- 
tional to T discussed in the previous subsection llll Al still 
prevails; (ii) for V >~ T the rising rate increases roughly 
linearly with V. 

On a finer scale, one may see small oscillations about 
the final value, especially in the curves for larger V. 
These result from the fact»i2i2& that in the presence of 

V the Kondo peak in steady state is split into two peaks 
at ±V/2 respectively. These split Kondo peak (SKP) 
oscillations were clearly identified in earlier workiiLi&S! 
They have a frequency almost precisely equal to V and 
a lifetime that is much longer than that implied by the 

V and T rate scales, something suggested by the pre- 
scient fact that in mean field theory the oscillations are 
undamped^ We put off a detailed analysis to a later 
section. 



C. Initial oscillations 

For small times, Yt < 10 in Figs. ^ and [3] there ap- 
pear oscillations, which are invisible on the scale of those 
figures. As opposed to the features mentioned in the 
previous two subsections, they are non-universal in the 
usual sense. However, for a given system, they take the 
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FIG. 3: Time dependent conductance vs. time. The ordinates 
are the time dependent conductance G(t) = I{t)/V in units 
of the open-channel conductance Go = 2(e 2 /2-kK). Here I(t) 
is the time dependent current, and V the bias. The abscissas 
are the time t after the gate switches the dot into the Kondo 
regime, in units of Both panels are for system one (see 

text), with T = 0.0015r (left panel) and T = 0.005r (right 
panel). The numbers in the legend are the biases V in units 
of r for the corrsponding curves in each panel. 



same form and frequency, independently of T and V, at 
least for values of these parameters that are much smaller 
than edot or T. They are shown in the top curve in the 
top panel of Fig. 0] for system one at T — 0.0015F and 
V = O.Oir. The next several curves show the initial oscil- 
lations for successively larger values of |edot|- The second 
panel more clearly establishes that the frequency of these 
oscillations is the difference between the Fermi level and 
the virtual level position. 




FIG. 4: Initial oscillations at short times. The top (long- 
dashed) curve in the top panel is the short time version of 
the system one (|eaot| — 2r) curve for bias V = O.Oir. Thus 
it represents the very short time behavior of the curves in 
Fig. The behavior is similar for all biases studied. The 
other three curves are for other values of |edot| whose ratio to 
F is indicated in the legend. In each case the time coordinate 
is rescaled as indicated, to show that the frequency of the os- 
cillation is about equal to the magnitude of virtual resonant 
level's energy (measured from the Fermi level). The renor- 
malized value of the latter is a bit larger than the bare value 
e c iot|. The bottom panel, showing the second derivatives of 
the curves in the top panel, more clearly identifies the period 
of oscillation. 



exact if the the current were instantly returned to zero 
at the end of the pulse. Since the switching goes to a 
non-Kondo region, the turnoff rate will be ~ T. Features 
that occur on that scale like the initial oscillations will 
be masked, but the V and T rates, as well as the SKP 
oscillations will be preserved. We have verified this by 
explicit calculations. 



D. Alternative roughly equivalent measurement 

The investigations described above shows that by mea- 
suring the transient currents in a quantum dot subject to 
a finite bias as a function of bias V and temperature, one 
can probe the bias-induced and thermal broadening of 
the non equilibrium Kondo problem. The most straight- 
forward implementation of such a measurement would 
be to measure the total charge transported through the 
dot Q(t) when the dot level is subject to a pulse train 
in which the dot level is raised to the Kondo regime for 
a time r. By measuring the total charge Q(t) trans- 
ported trough the dot during a single pulse as a function 
of pulse length t and taking the derivative with respect 
to r, dC i , a quantity is obtained that closely follows the 
transient currents in the dot. The equivalence would be 



IV. ANALYSIS AND INTERPRETATION OF 
ZERO BIAS RESULTS 

A. Large time limit 

In Fig. we display the large time limit of our calcula- 
tions for systems one and two. Since the bias is essentially 
zero, these represent the steady state linear response con- 
ductance. This figure shows that our results agree with 
general expectations including approximate universality. 

The exact asymptotic curve at large ln(T/Tk) 

2. = —^. (4) 

G 161n 2 (T/T K ) 

was first calculated by Abrikoso\i2£ for the Kondo impu- 
rity problem, and has been more recently adapted and 



5 



applied applied to quantum dotsS It can also be de- 
rived by the so-called poor man's scaling methodi^Si^i 
However, if one applies the perturbative procedure used 
in Ref. 21 for large V to the similar large T case, one 
finds that within NCA, the asymptote is 4/3 times the 
value in Eq. The NCA asymptotic is the one shown 
in the figure. Our lowest temperature point (not shown) 
is slightly above the unitarity limit, another NCA error 
also found by othersiSI 
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FIG. 5: The large-time limit of G(t)/Go vs. temperature. 
Since V ~ 0, this is the equilibrium value of G/Gq. The 
points labeled SI and S2 are from system one and two, respec- 
tively (see text). The curve labeled asymptote in the legend 
represents the large T asymptote within NCA (see text). 



B. Extraction of the rise rate 

Here we use previous analytic results^ to extract the 
rise time r from the data of Fig. 2] What was found there, 
was that the current response to a stepped turning on of 
the Kondo interaction would be the same as the equi- 
librium response to a time-dependent effective tempera- 
ture T c given by T c = T cothirTt/2h. We start with the 
tautology SI /I = [d(ln G)/d{\nT c )} 6T C /T C . Defining the 
fraction / as the finite difference 



and 



/ 



fo 



J(o°) - Ht) 
J-(oo) : 



we have in the large t limit / 
which becomes 



d(lnG) 
d(lnT)' 

fo(cothirTt/2h 



I 



/oexp ( -- 



(5) 

(6) 
l)/2. 

(7) 



= ttT 



(8) 



This suggests that we should fit the upper part of the 
our curves in Fig. ^ to Eq. I|7|. and compare the resulting 
r with Eq. (JSJ. Of course, the derivation in Ref. 1 1(3 was 
only to lowest logarithmic order, and certainly cannot be 
expected to be valid for T <^ Tk, and we might also 
expect the possibility of logarithms of T/Tk to become 
predominant at large T. Nevertheless, we find no evi- 
dence for the latter effect in the curves of Fig. ^ We 
display the rates (= 1/r) that we find in Fig. HJ1 The 
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FIG. 6: Rate (1/r) at which the conductance approaches its 
final value versus temperature T, with zero bias V = 0. The 
solid dots are for system one (see text), while the open circles 
are for system two, which has a different Kondo temperature 
Tk. The solid curve is a straight line of slope n through the 
origin (see Eq. ||HJ). 

quite good agreement of the curves in Fig. Hoover a major 
part of the range is satisfying. The curves also strongly 
suggest that the low temperature limit of Eq. © replaces 
T by something of the order of Tk- Our points appear to 
deviate from Eq. (JSJ at higher temperatures. However, 
our fitting procedure is less accurate at high tempera- 
tures where G(oo) is small, and our estimated numerical 
error is of the same magnitude as this deviation. For this 
reason, we are not prepared to say with certainty whether 
this is a real effect or not. An analysis of the temporal 
evolution of the instantaneous spectral function of the 
dot level does reveal significant spectral reshaping of the 
Kondo resonance well beyond the time where the conduc- 
tance has saturated^ This effect suggests the presence 
of a slower time scale at large temperatures. We might 
certainly expect a slower timescale to emerge, if we apply 
the arguments of Ref. to the low bias, high tempera- 
ture case, as done in Appendix ^ A more detailed in- 
vestigation of the origin of this effect and the possibility 
for its experimental detection is in progress. 



6 



V. ANALYSIS AND INTERPRETATION OF 
FINITE BIAS RESULTS 

A. Large time limit 

We begin by plotting the large time values of G(t) in 
our finite V calculations. This of course is the steady 
state conductance, and our curves in Fig. {7\ supplement 
results previously calculated^! for T = 0, which are in- 




FIG. 7: Long time or dc conductance vs. bias. The ordinate 
G/Go is the dc conductance in units of the open channel con- 
ductance Go — 2e 2 /2-kTi for the two spins, while the abscissa 
is the bias V in units of the Kondo temperature Tk- The tri- 
angles represent the data of Ref. |^ with the abcissa rescaled 
to correct for the ~ 30% difference between the two defini- 
tions of the Kondo temperature. The open circles, squares, 
and diamonds represent the long time saturation of our time- 
dependent conductance curves for system one at the temper- 
atures indicated. The smaller solid symbols represent the 
values at the same respective temperatures for system two, 
which as a different Kondo temperature (see text). 

eluded in the plot. Since the formula for the Kondo tem- 
perature used there gave a value some 30% greater than 
Eq. (J2J, their values for V/Tk were rescaled appropri- 
ately, so they could be put on the same graph with ours. 
Our values show approximate universality as discussed 
earlier and agreement with Ref. in the appropriate 
regions. 



B. Extraction of the rise rate 

The very large time behavior of our curves in Fig. 
is dominated by the decay of the SKP oscillations, ft 
is clear that their amplitude is very small and for the 
smaller V^'s entirely negligible, and they will have no 
effect on any common sense definition or experimental 
measurement of rise time. Despite the SKP oscillations, 



the form Eq. J7J) well describes the upper part of the time- 
dependent conductance curves in Fig. [3] up until they 
reach within ~ 1% of the saturation value. Therefore, 
as a practical empirical technique to characterize the rise 
time, we use the same technique as for V = 0. The results 
of this analysis are shown in Fig. |H1 




V/T K 

FIG. 8: Rate (inverse time) at which the conductance ap- 
proaches its 'final value' versus bias V, at several tempera- 
tures as indicated. The open circles, squares, and diamonds 
are for system one (see text), while the corresponding solid 
symbols are for system 2, which has a different Kondo tem- 
perature. Our V = points are shown in the figure at 
V/Tk — 0.1. The solid lines are the predictions of Eq. JjJJ. 



C. Interpretation 

If the source of the T in Eq. © is the same as 
the source of T in the Korringa rate, Eq. IjAljl . then 
it arises because a change in S involves the absorp- 
tion of an electron-hole pair of zero energy from the 
leads. The phase space for such a process is doubly re- 
stricted by the Pauli principle, and produces the factor 
J def(e)(l — /(e)) = T where /(e) is the Fermi function. 

The application of these ideas to a biased dot allows 
the phase space factor to be tuned in a continuous fashion 
by varying the bias V across the two leads of the dot, 
which for simplicity we assume to be symmetric under 
lead interchange. The phase space restriction factor in 
this case is given by \ Y,i,i< I de //( e )(! ~ Z'( e ))> where 
the indices I and V designate which lead is referred to. 
For example, the Fermi functions /i(e) and /2(e) have 
Fermi levels displaced by ±eV/2, respectively, where e 
is the magnitude of the electronic charge and V is the 
bias. The integral above can be evaluated analytically, 
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with the result that Eq. (JSJ) should be replaced by 



where 



F(*) = i ll + l + r 



xe 



1 



(9) 



(10) 



In writing Eq. (|l(Jf) we have sacrificed conciseness to fa- 
cilitate clarification of the origin of the terms in paren- 
theses, in terms of the wave function of the annihilated 
particle-hole pair. The first two terms arise when the 
two components of this object are in the same lead; in 
this case the existence of V has no effect on the result; 
the phase space for these processes is still constricted, 
and the contribution to 1/r is still small. For the third 
term, the particle is on lead 1 and the hole on lead 2; 
here the phase space is opened wide by V. Finally for 
the fourth term, the particle is on lead 2 and the hole 
on lead 1; the phase space is, aside from an exponen- 
tially small tail, closed off entirely as V is increased. So 
the essential physical feature deriving from Eq. (|10|) is 
that the factor of T in Eq. JSJ is replaced at large V by 
jV. The notion of an expanded phase space is implicit 
in the Anderson model rate calculation of Wingreen and 
Meir— and in a different context in the work of Kaminski 
et alw^k In any case, the rate of Eq. provides a way to 
rationalize the calculated points in Fig. [8] in a parameter 
free way, which has predictive power for the rise-time of 
the conductance. The comparison shown in that figure 
shows that it captures the main trends of the NCA re- 
sults, although not with such good agreement as for the 
V = case (Fig.©. 




FIG. 9: Time-dependent spectral functions for the time scale 
relevant for the rise time of the conductance. These curves are 
for system one with T = 0.0015r = 0.69T K , and V = 0.08T = 
37Tk. The curves are at equal time intervals At — 2.8r~ 1 = 
0.006T" 1 = 0.17(4/7T V), where the 4/ttV factor is motivated 
by the high temperature asymptote of Eq. @. This factor 
has the value of unity at t = 6 At (thick solid black curve). 
The lowest curve is for time (t = 2 At) after the virtual level 
parameter was switched to it's final value of — 2F. By curve 
15 the conductance has for all practical purposes reached its 
dc value. 



probably unavailable through perturbation theory from 
the true steady state, and hence will not likely appear in 
steady-state correlation functions. 



VI. TIME-DEPENDENT SPECTRAL 
FUNCTION— TWO TIMESCALES 



B. Slower timescale 



A. Faster timescale 

As an aid to the interpretation of the conductance 
behavior, we display the time-dependent spectral func- 
tions A<iot(v,t) for the dotS Fig. shows the case 
T = 0.0015r = 0.69T K and V = 0.08T = 37T K . The 
snapshots shown here are for the time scale appropriate 
for the rise time of of the conductance, as shown in detail 
in Fig. ^3 Fig. [5] shows the rise of the Kondo peak as 
a mostly smooth structure of half-width ~ V, with the 
individual split peaks mostly undeveloped. This struc- 
ture appears to be converging toward a quasi-stationary 
steady state. It is clear that this timescale, intermedi- 
ate between the trivial very rapid timescale governed by 
the T rate, and the longest timescale yet to emerge, is 
the one governing the rise time of the conductance. It is 
probably only available through a fully non-equilibrium 
theory, since it is the scale relevant for approaching a 
steady state that is metastable at best. It is therefore 



Fig.llllshows the development of the spectral function 
on the longest timescale, showing the rate at which the 
individual split peaks emerge from the quasi structure of 
width V formed at the metastable point. For V » Ik, 
this latter rate is much smaller than that illustrated in 
Fig- El i n this case by an order of magnitude. During the 
whole course shown in Fig. 1111 the conductance changes 
only by the small oscillatory amounts indicated by the 
SKP oscillations. Indeed, the area under ^4dot (when 
T <C V) in between —V/2 and V/2 is a virtual con- 
stant in time in this time regime, with the additional 
area under the split peaks being almost exactly compen- 
sated by the loss in area between peaks, continuously in 
time. Indeed, we have verified that the conductance in 
this regime, which still has a tiny fluctuation due to the 
SKP oscillations, is given by the steady-state formula (for 
example, Eq. (12) of Ref. ll9l) . provided that the ^4dot used 
is the time-dependent instantaneous spectral function. 32 . 
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FIG. 10: The rise in conductance as the Kondo peak forms for 
bias V = 0.08F at T = 0.0015. The solid curve is a magnified 
section of the corresponding curve in Fig.[3]with the horizontal 
axis shifted as indicated in the legend. The dashed curve is 
the indicated exponential fit. The rate R obtained in this way 
characterizes the rise time of the conductance, and this is the 
rate used in the preparation of Fig. |H] Of all the curves in 
Fig. [3 this is the worst case for the exponential fit, because 
here the SKP oscillations are largest. The first period of these 
oscillations is partially visible here. Multiple cycles are shown 
in the larger magnification of Fig. 1121 



VII. SKP OSCILLATIONS AND THE LONGEST 
TIME-SCALE 

A. SKP oscillations in the quantum dot 

In Fig. ^3 which extends a little beyond the range of 
Fig- El one can see the beginnings of the SKP oscillations. 
The continuations of them on a highly magnified scale 
are shown in Fig. 1121 Generally, many complete periods 
of oscillation can be seen, and this period corresponds 
almost exactly to an angular frequency of V , becoming 
closer and closer from below as V gets larger and larger. 
The details of this analysis is given in Appendix^] There 
it is also shown that the decay rate is close to constant 
(exponential decay) . For large V this rate is substantially 
smaller than the rates identified earlier characterizing the 
rise time of the conductance. The method of extraction 
of the rates is described in Appendix [BJ 

For the case shown here, this rate is on the order of the 
distance to the inflection point on either of the peaks in 
the converged spectral function^ The rate is not incon- 
sistent with, but possibly a little slower than, the slow- 
scale rise rate of the time dependent spectral function, 
although the latter is more difficult to pinpoint. We can 
with much less ambi guit y compare this rate with the 27 
rate identified in Ref.|2l], and calculated for T = 0. What 
we find is that if we divide the SKP oscillation decay 



FIG. 11: Time-dependent spectral functions for the time 
scale relevant for the development of the split Kondo peak 
and the damping of SKP oscillations. This is a continuation 
of the snapshots shown in Fig. |U] so all the parameters have 
the values shown in the caption of that figure. Here, how- 
ever, the time interval between the successive snapshots is 
much greater, varying from 5 At to 50 At, as indicated by the 
labeling in the legend. 

rates by two, then these values at our lowest tempera- 
ture (~ 0.3Xk) agree within their accuracy with the 27 
curve. This is shown in Fig. 1131 We are not certain why 
half the decay rate is what seems to correspond to 27, but 
presumably the issue is dephasing which occurs both at 
—V/2 and +K/2 giving roughly additive contributions. 
Another possibility is that T = (which is currently un- 
available for us) will be further from our lowest points at 
T = 0.34Tk than we expect. But a definitive answer will 
have to await further and more complete studies. What 
is very clear, however, is that the damping of SKP oscil- 
lations is controlled by a rate much slower than ~ V for 
large V, and which is numerically quite close to the 27 
rate. 



B. Decay of SKP oscillations in A- fold degenerate 
models 

We conclude this section with a connection with mod- 
els which increase TV from the N = 2 appropriate for 
quantum dots. Repeating the large V NCA perturbative 
analysis 21 for that case gives 



As noted in Appendix the NCA does not need cor- 
recting by multiplying this by S(S + 1), which would be 
proportional to N 2 at large N. Eq. (|llfl indicates that 
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FIG. 12: SKP oscillations for bias V = 0.08r at T = 0.0015. 
The solid curve is a magnified section of the corresponding 
curve in Fig. [3]with the position of the horizontal axis shifted 
as indicated in the legend. The dashed curve is a progressivly 
more magnified version of the solid curve, so that more SKP 
oscillations can be seen. As indicated in Fig. 1141 the angular 
frequency of these oscillations is almost exactly equal to V. 
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FIG. 13: Decay of SKP oscillations vs. bias V. The empty 
squares, diamonds, and triangles represent data derived from 
the damping rate of the SKP oscillations in our calculations 
of conductance vs. time at the temperature indicated respec- 
tively. The solid circles are the calculations of 27 at zero 
temperature in Ref. 12 lL 



the decay rate of SKP oscillations approaches zero as 
1/N for large N. This is consistent with the trend of 
NCA studies^ of this damping as a function of N, and 
with the faciei that these oscillations are undamped in 
mean field theory (N — 00). In this case the effect of the 
vertex correction (see Appendix EJ is dramatic indeed, 
making a qualitative change, as opposed to the ~ 30% 
correction for the N = 2 quantum dot case. 



VIII. ADDITIONAL CONCLUSIONS 



Many of our important conclusions arc contained in 
the abstract and final paragraph of the introduction. We 
mention here the few that are not. 

First, for V ~ the rise-time of the conductance is 
very accurately characterized by the simple expression 
7rT, Eq. JHJ, over a wide range of temperature T > Tk, 
and appears to heal towards Tk for smaller T. Second, 
for V > 4T, one should replace T by V/4 in the above. 
These rates appear to lead toward a quasi-metastable 
point, and do not negate the emergence of slower rates 
at longer times, and which in the case of larger V con- 
trol the damping of SKP oscillations. Third, at the very 
beginning of switching the gate, there are small initial 
nonuniversal oscillations at a frequency corresponding to 
the dot level's separation from the Fermi level (for V's 
small with respect to that separation). Finally we pre- 
dict for large V ^ Tk for A^-fold degenerate models that 
the damping of SKP oscillations decreases as 1/N, 
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APPENDIX A: THE KORRINGA RATE 

A half a century ago, the framework and starting point 
for discussion of the time scale r for a localized spin in 
an electron sea was set by Korringa, 35 whose contribution 
has been promulgated at the textbook level for decades^ 
In simple terms 1/r gives the fractional rate at which a 
component of spin S representing a magnetic impurity 
(or a quantum dot) is changing due to the electrons in 
the conduction band (or in the leads of a quantum dot). 
It is given by an expression of the type 



= aT, 

7"kor 



(Al) 



where a is a dimensionless constant, and T is the temper- 
ature. Although Korringa's original derivation applied 
to a nuclear spin where the interaction with electrons is 
dipolar, it has been widely applied to impurity electron 
spins as well, using the Kondo model 



(A2) 



for the interaction. In Eq. ljA2fl , ki is a generalized quan- 
tum number for the i th conduction (or lead) electron, 
is its energy, and is its spin operator. For a symmetric 
quantum dot, k is assumed to include the information 
on which lead is referred to. The exchange coupling J is 
taken to be independent of fc, aside from the usual high- 
energy cutoff at energies further than D above or below 
the Fermi level. In terms of Anderson model parameters 
J = — 2\V\ 2 /edot, where here V is the matrix element 
in Eq. QJ, an d not the bias voltage. To lowest order in 
J, the quantity a in Eq. (HTJ is given, for example, by 
Eq. (1.16) in Ref. OH 



a = n(JpY 



(A3) 



where p is the density of states in the leads (both leads 
together), when applied to a symmetric quantum dot. 

One should note that (|A3|) applies to a general spin (5) 
Kondo model, (5 — | for a quantum dot), and there is 
no factor of 5(5 + 1). For 5 = \ this can be shown to be 
incontrovertible by the solution of a simple master equa- 
tion for the occupation probability of up or down spin 
states P a = —R(P a — P- a ), where R = R<^ = R^ is the 
spin-flip rate of the dot spin due to it's interaction with 
the lead electrons. This is readily solved for the decay 
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rate of an averaged spin component (S m ) = — (S' m )/rkor 
where l/rk r — 2i? = aT, where was evaluated by 
the Fermi golden rule, giving (|A1|I with a given by fO . 

While it is true that the equilibrium pseudo-fermion 
propagator's self-energy's imaginary part does indeed 
contain the S(S + 1) factor [Eq. (C4) in Ref. the 
time decay rate of the spin correlator does not, because 
of vertex corrections, which are summed in Ref. |33 via 
the Kadanoff-Baym equations. 

When J is rescaled via the poor man's scaling tech- 
nique, Jp —>■ — l/ln(T/XK), as done by many authors to 
obtain results valid for ln(T/Xk) ^> 1, one gets 



1 



ttT 



In 



2 T 



(A4) 



ima for the case shown in Fig. It is obvious that a 
fixed period is maintained over many cycles. This period 
is determined by a least squares fit to the minimum and 
maximum positions. The zero crossings straddle the fit- 
ted line because the peak-to-peak envelope center is still 
a small way below its final value. The decay rate was 
determined as shown in Fig. 1151 by a least square fit to 



Similarly one gets large V results through the replace- 
ment T — ► jV in Eq. I|A4|I . These agree exactly with 
what one gets from high V or high T expansions^ of the 
NCA equations, and we conclude that NCA gets the right 
answer here. Of course, the NCA misses the factor | in 
the conductance, Eq. (@J, as previously pointed outiSi 



APPENDIX B: DETAILS OF SKP OSCILLATION 
PARAMETERS 
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In Fig. E| we map the successive positions of the zero 
crossings as well as the positions of the maxima an min- 
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FIG. 14: Positions of extrema and zeros of SKP oscillations 
for V = 0.080r as in Fig.^U The slope of the fitted line gives 
a frequency of 0.0783P. The frequency is always a speck less 
than V, but for our largest V's, we find the difference to be 
only around 0.1%. 
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FIG. 15: Determination of the decay rate of the SKP oscilla- 
tions from the positions of the minima in Fig. 1121 The slope of 
the curve gives a rate of 0.015P for V = 0.08r and T = 0.0015. 
These rates are universal functions of V/Tk and T/Tk, and 
hence depend on T, E, and D only through Eq. In the 
case of this figure, T K = 0.0022r. 



the logarithm of the difference between the dc conduc- 
tance and its value at the minimum points vs. the time 
corresponding to these points. If one looks closely, there 
appears to be a small curvature at the beginning. This is 
probably real and not noise, as according to Fig, im the 
split Kondo peaks in the instantaneous spectral functions 
are still narrowing slightly in this region. Aside from re- 
jecting early points that were obviously out of line, we 
did not attempt to account for this in the analysis. The 
small misalignment of the final point, on the other hand 
is certainly due to noise, an ubiquitous feature for ex- 
tremely large times, which deterred us from trying to 
distinguish between the 'initial' and 'final' slopes. 



